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The dynamics of the gauge vortex with arbitrary form of a contour is considered in the framework 
of the nonrelativistic Abelian Higgs model, including the possibility of the gauge field interaction 
with the fermion asymmetric background. The equations for the time derivatives of the curvature 
and the torsion of the vortex contour generalizing the Betchov-Da Rios equations in hydrodynamics, 
are obtained. They are applied to study the conservation of helicity of the gauge field forming the 
vortex, twist, and writhe numbers of the vortex contour. It is shown that the conservation of helicity 
is broken when both terms in the equation of the vortex motion are present, hrst due to the exchange 
of excitations of the phase and modulus of the scalar field and the second one due to the coupling 
of the gauge field forming the vortex, with the fermion asymmetric background. 
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1. Introduction. The string-like (or vortex) solutions of the field-theoretical models are widely discussed in 
application to numerous physical systems, from the really observed in condensed matter physics (quantized vortices in 
He'^ [l|, Abrikosov lines in type II superconductors 0, vortices in Bose - Einstein condensates Q etc.) to hypothetical 
cosmic strings [J . Of particular interest are the situations when the dynamics of the vortex contours can be deduced 
from the field equations of the underlying field theory. Different type of models are invoked to study the dynamics. In 
particular, the dynamics of cosmic strings was studied in the framework of the Nambu-Goto action ii. The latter 
was shown to result Q from the relativistic version of the Abelian Hires model (AHM). The solution in the form of 
the static cosmic string was obtained from the static version of AHM Q . Another limiting case is the nonrelativistic 
one, which can be deduced from the nonrelativistic Abelian Higgs model (NRAHM). The static limit of this model 
coincides with the Ginzburg-Landau theory Q . 

The vortex in the models with the local gauge symmetry is called the local one. The equation of the nonrelativistic 
motion of the local vortex, in the idealized situation of the negligible dissipation, was obtained in the work llOll in the 
framework of NRAHM. The feature of the equation of the vortex motion in the model considered in Ref. is that 
except for the usual term stating that the vortex ring velocity points to the direction of bi-normal, there appears the 
additional term arising due to the exchange of the excitations of the phase and modulus of the scalar field between 
different segments of the contour. Yet the model considered in is not the most general one. As was pointed out 
in m , there could be another term in the gau ge v ortex equation of motion in case when the gauge field is coupled to 
the asymmetric background of chiral fermions |I2l-lbM. In the static case, the influence of this additional term on the 
form of the vortex contour was considered in Ref. [ll|. 

Another interesting aspect of the time-dependent gauge vortex solutions is that they can elucidate the interrelations 
among such characteristics of the closed contours as the helicityof gauge field forming the vortex, the writhe and torsion 
(twist) numbers and their possible dependence on time |18l - l20l |. In case of pure hydrodynamics, the corresponding 
equation was deduced in Ref. [2l|. Note that the problem of dynamical evolution of the vortex filaments has a long 
story traced to the beginning of the twentieth century . The corresponding equations for the time derivatives of the 
curvature and torsion were rediscovered, in particular, in Refs. [221 l23l|. The historical review of these developments 
is presented in Ref. [^ . 

The aim of the present work is to consider the dynamical evolution of the gauge vortex string in NRAHM, including 
the possibility of its interaction with the static fermion asymmetric background. 

2. Nonrelativistic Abelian Higgs Model with the gauge vortex. Nonrelativistic Abelian Higgs model 
incorporating gauge vortices is given by the following Lagrangian density 
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where E = —dtAjc—'^tf and H = 'V x A are the electric and magnetic field strengths, the four-vector gauge 
potential is = {(p,A); q, m are the charge, mass, of the particles forming the condensate of the scalar field, c is 
the velocity of light. The quantities uq and po are, respectively, the density of the scalar field condensate and the 
homogenous positive charge density introduced to provide the net neutrality of the system: 


Po + qno = 0. 


( 2 ) 


The coupling constant g can be related to the sound velocity Cs- See Eq. © below. The limitire case of this model, 


with the gauge coupling constant set to zero, corresponds to the model of Gross-Pitaevskii [2^, [2^ possessing the 


solutions in the form of global vortices observed in Bose systems like Hell or BEG. The vortex is represented as the 
line of singular phase Xs of the scalar field 'ip = which forms the spatially homogeneous condensate of the 

density uq everywhere except for the vortex core, where it goes to zero at the transverse distances of the order of the 
healing (or correlation) length 

The notations in Eq. o are as follows. The four-vector 


a-tj. = - dfj,Xs 

q 

is the four-gradient of the singular phase, 

[V X V]xs = J - X), (3) 

serving as the source of the gauge vortex with the unit flux quantum whose space-time location is given by the vector 
X = X{t,a). The length of the contour, I = \X'\da, is the natural choice of a. The explicit expressions for 
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found from Eq. ©, in the gauge Va = 0, are: 
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Hereafter, the prime (dot) means the differentiation with respect to the contour parameter a (time), and [a x b] 
[(a, 6 )] stands for the vector (scalar) product of two vectors a and b. 

The term oc 5 in Eq. © can be rewritten in the form 
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demonstrating that the first term in the right hand side describes the s-wave interaction between the bosons forming 
the condensate, the second one is with the chemical potential /r = griQ pertinent for the Bose gas with the s-wave 

interaction [^ . and the third one is irrelevant constant. After integration over space and going to the Hamiltonian 
this second term goes to —gN thus showing that the actual Hamiltonian is in fact H' = H — gN This chemical 
potential /i is due to the background of the bosonic condensate. The equation of motion of the gauge vortex resulting 
from NRAHM was shown [l^ to look like 


where 
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is the the London penetration depth, 
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is the velocity of sound [^ . 




h 

2mcs 


( 7 ) 

( 8 ) 


is the healing length (the characteristic size of the vortex core). The term oc l/Cg and the intermediate scale 


As = Al — <c Al, (9) 

c 

as was shown in Ref. are generated dynamically, due to the exchange of the fluctuations of the phase and the 
modulus of the scalar field ip between the different segments of the gauge vortex. The term [X x X'] originates 
from the variation of the term J q\ip\‘^aodtd^x in the action 0 - Eor further use let us introduce also the quantum of 
magnetic flux 


^ 2Trhc 
$0 =- 

g 


( 10 ) 


It is known that the presence of the background of chiral fermions (whose asymmetry is characterized 

by the chemical potential /r/) coupled to the gauge field (in particular, the Abelian one [laB), induces the term 
—gfNcs in the Hamiltonian, where 


^CS = ^ y d^x{A, [V X A]) (11) 

is the Chern-Simons number carried by the Abelian gauge field of the vortex. Coupling of the Abelian gauge field 
with the chiral fermion [l5l| takes place, in particular, in the standard model of particle physics 0 where the right 
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electron is coupled to the hypermagnetic gauge field. As is seen from Eq. (ED, the Chern-Simons number in the 
Abelian case is proportional to helicity, 


hA= J d^x{A, [V X A]). 

Using the expressions for A and B = [V x A] from Ref. [la, 
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one obtains the expression for the helicity of the single vortex with the unit flux quantum and its variation over X: 
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where Xi ^2 = -X'(tTi^ 2 ), X 12 = Xi — X 2 . Since the segments of the vortex contour located at the distances greater 
than A l give exponentially damped contribution, one can use the expansion in z = cr 2 — ci, 


2 3 

X[a + z)= X{a) + zX'{a) + ^X"(ct) + ^X"'(a) + ... 

2 D 

and keep only the first non-vanishing term surviving in the limit Al —>■ 00 : 

( 1)2 r (f ,2 f 
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The variation of the corresponding term in the action, AS = J dtNcs , gives the term oc [X' x X'"] in the equation 
of motion m which now looks like 


[X X X'] = 7 X" + To 
where the following shorthand notations are used: 


[X X X'] X X' 
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The term 7 X" is due to the tension which push the vortex segment to the direction of the normal vector n. Other 
terms are discussed below. When both To and /I are set to zero, the known relation X = 'ynb is recovered. In this 
case the static solution with zero velocity is possible for the straight vortex for which k = 0, t = 0. The particular 
dynamical solution to the latter equation describes the vortex ring moving with the velocity V = •jk in the direction 
orthogonal to the ring plane. If Tg = 0 , /i 0 then the nontrivial static solution in the form of the helix with k =const 
and T = "f /fl\s possible [HI- 

Taking into account the fact that the physical motion of the vortex is orthogonal to the unit vector X', one can 
rewrite Eq. ED in the form of the wave equation with the additional terms: 


4-X -X" + -[X X X'] T^X' X X'"] = 0, 
Co 7 7 
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where 
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The interpretation of the additional terms is the following. The term oc [X x X'] is the analog of the Magnus force 
acting on the moving vortex due to the nonzero circulation of the supercurrent. In the classical hydrodynamics, in 
the direction orthogonal to the displacement velocity, the velocity of the circulating flow adds with the displacement 
velocity on one side, and subtracts from it on the opposite side. By Bernoulli’s equation, the pressure in the liquid is 
greater in the regions where the velocity is lower, yielding the force in the direction orthogonal to the vectors X and 
X'. Of course, in the case of our interest the situation is much more intricate because one should take into account 
the role of a number of quantum excitations, but the qualitative picture is probably the same. 

The term oc [X' x X'"] is induced by the anomaly. Indeed, the variation of the induced Chern-Simons term in 
the action over vector potential A results in the anomalous current jan oc /r/iT pointed along the direction of the 
magnetic field. This relation is the corner stone of the widely discussed chiral magnetic effect. See the review [l^ and 
references therein. This current is subjected to the action of the Lorentz force oc [jan x H] from the magnetic field 
produced by the nearby segments of the gauge vortex. Taking into account the exponential damping of the magnetic 
field strength and the expansion m, one can write the leading contribution to the anomalous part of the Lorentz 
force per unit length: 


/ OO 

dz[X'{a) X X'{a + fif[X' x X"']. 

-OO 


( 21 ) 


So, the discussed term is the macroscopic manifestation of the quantum anomaly. The helical form of the vortex 
contour in the static case is due to the balance between the tension and the anomaly-induced Lorentz force. 

3. Time derivatives of the basic contour characteristics. As is known, the derivatives over the contour 
parameter a of the basic unit contour vectors of normal n, bi-normal b, and tangent X' constituting the right triple, 
X' = [n X 5], are given by the Frenet-Serret equations 
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where k, t are the curvature and torsion, respectively. One can write down the analogous set of equations for the time 
derivatives of the above vectors taking into account the fact that they are unit ones and that the physical motions of 
the contour are transverse (i.e. orthogonal to X'): 
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Equating the commutator of the derivatives over contour parameter and time to zero one finds that the matrix A in 
(1^51) is antisymmetric, abn = —dnb, a,nx = —dxn, and the following equations are valid: 

^nb — T ^Xb}/ 

k = — Taxbj 

f = naxb + aJib, (24) 

so that 

( 0 axn axb \ 

^Xn 0 ^nb j (25) 

— axb —anb 0 J 


One can relate k and r with the components of the vortex velocity X, or, equivalently, with the dual vector W = 
[X X X'] = WnTi + Wbb, X = Wnb — WbTi. Differentiating the definition of the vector W over a and using the 
Frenet-Serret equations and Eq. (1^ . one obtains 
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This equation permits one to make the identification 

axn = rWn, 

axb = < - rWb, (27) 

so that 


k = -{Wl + TWr^y - - TWb), 

T = k{W^ - TWb) + {[« - TWb)' - t{WI, + TWn)] /«}' ■ (28) 

With these ingredients one can rewrite Eq. (1171) in terms of the components of W : 

Wn = k(7 - /It) + To {Wb + [« - rWb)' - t/W/ + rW„)] W„/k} , 

Wb = fin’ - To - [« - TWb)' - t{WI + TWn)] Wb/n} . (29) 

The system of equations (l28l) is in fact the result of definitions which are not related to the specific dynamical model. 
The latter is provided by Eq. (17^ . In particular, when Tq and fi both vanish, Wn = 7 k, Wb = 0, one finds in this 
limiting case that 


k = — 

T = 7 


^{ 2 kt + kt'), 

/ f 1^" 

KK +- 



(30) 


After the re-scaling k —>• t —>• 7“^/^r, a —>■ 7^/^cr the parameter 7 drops from equations, and one gets the 

Betchov-Da Rios equations (2ll - l^ known in hydrodynamics for a long time. 

Equations (I29p are nonlinear, so let us find their solution which is exact in the fermion chemical potential fl but 
perturbative to the first order in Tq. The zeroth order solution is 

Wl;°'> =jiK'. (31) 

The corresponding analog of the Betchov-Da Rios equations found from (1751) looks as 
= —(7 — /ir)(2K'T -I- kt') — JI/k" — kt'^), 

= kIjk' — fl{2K'T + kt')] -(- { [(7 — jlT){K" — KT^) -|- P,{ 2 k't + Kt'Y] /k} (32) 

It is straightforward to obtain the first order correction oc Tq using the expressions 
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, = -To <{ - - [(W;, - TWb)' - r{Wl + TWn)] Wb 

Hi 
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(33) 


but the corresponding expressions are excessively long. Instead, one can find the time derivative of the total torsion. 
From the third equation in (IMl) and the second equation in (1771) one obtains 


ida = Kaxbdcr = — ^(k'IE„ + KTWb)da. 


(34) 


Inserting here Eqs. (15T1) and (155)) and using Eq. (1571) . after some lengthy algebraic manipulations and integrations by 
parts, one finds 


/ 


ida = - 7 T 0/2 (p K T da. 


(35) 


This expression shows that the total torsion of the gauge vortex in the nonrelativistic Abelian Higgs model is not 
conserved in case of the presence of coupling of the gauge field forming the vortex, with the fermion asymmetric 
background described by /i 7 ^ 0. Coming back to the physical quantities one obtains 
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One should have in mind that the equation of motion ([5]) was obtained in the paper 0 in the London limit ^ 1, 
Ag/^ » 1, so Eq. (15^ tells us that the time derivative of the total torsion, being proportional to the small factor l/c^, 
is small but not negligible, because it is multiplied by two large logarithmic factors. Note also that the coherence 
length ^ characterizes the radius of the vortex core where the condensate of the scalar held vanishes. 

The time derivative of helicity is obtained from Eq. (HU): 
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where Tw is the twist number. 
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Tw — '''d.a = — j){X', [n' x n])da. (38) 

One can relate the above quantities with the writhe number. The term ([AC' x X"'],5X) in the integrand of 
expression for the variation of helicity (IT6l) . hence the corresponding term ([X' x AC], X'") in the expressions for the 
time derivatives dSH), has the geometrical meaning. Indeed, let us hnd the writhe and its variation: 
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The ambiguity 0 x 00 in the integrand can be regularized with the help of expression 

-x^/2e 
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5^0 (27r)3/2^3 

Then, using the expansion (1151) in z = a 2 — cri, the regularized expression for the variation of the writhe number can 
be written in the form: 
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Combining this equation and Eq. dSZl) one can see that 


(41) 


-^(Wr + Tw)=0, (42) 

at 

i.e., the sum of the twist and writhe is conserved as it should (Tsl - f^ 1^. In the meantime, the helicity of the 
gauge field configuration in the model considered here is not conserved. It is essential that the conservation of helicity 
is broken due to the fact that both terms in the equation of the vortex motion are present, first due to the exchange 
of excitations of the phase and modulus of the scalar field resulting in the factor oc c“^, and the second one due to 














the coupling of the gauge field forming the vortex, with the fermion asymmetric background oc . The latter serves 
as the source and sink of the helicity. Such coupling is in fact the consequence of the U(l) anomaly of the current of 
the chiral fermions [H. This anomaly is a quantum phenomenon hence the appearance of the Planck constant in the 
right hand side of Eq. (1551) . 

4. Discussion and Conclusion. The vortex solutions with nontrivial dynamics in the gauge field models like the 
nonrelativistic Abelian one considered here, are interesting in that they allow to apply some general mathematical 
concepts like the topology of curves and the equations for the time derivatives of the curvature and torsion, to the 
situations beyond the earlier considered hydro- or electrodynamical ones (Tsl - f^ . As is shown in the present work, 
one can relate the above time derivatives with the vortex velocity. The latter, in turn, is governed by the equation 
of motion specific for a given underlying model. In the presence of both the dynamical correction to the equation 
of motion due to exchange of the fluctuations of the phase and modulus of the scalar field, and the coupling with 
the fermion asymmetric background, the helicity of the gauge field is not conserved. Contrary to the situation in 
magnetohydrodynamics where the helicity conservation is violated by the dissipative effects (finite conductivity of the 
medium), in the present work, the nonconservation is reversible, because both signs in the right hand side of Eq. (13611 
are possible. 

As for the application to physics, one can hardly hope to observe the closed or curved vortices in type II supercon¬ 
ductors immersed into the static external magnetic fields. However, their production could be quite possible in the 
fast temperature quench See, for example. Ref. [s^, where such a mechanism was considered in the framework 
of the 2D Abelian Higgs Model. Also, numerous simulations of the phase transitions in the early Universe point to 
a considerable portion of the closed cosmic strings j^. The type of the model considered in the present work is one 
where the charged particles interacting with the gauge field, condense uniformly except for the vortex core. The total 
neutrality is provided by some static background of the opposite charge. The cosmic strings of such type are not 
excluded in the models of the dark (hidden) sector [H, [s^. However, the relation of the present work parameters q, 
no, and g with ones from the models of the hidden sector, goes beyond the scope of this paper. 

The global vortices, i.e. the vortices in the models with the global U(l) symmetry, of arbitrary form, were observed 
in the dilute Bose or/and Fermi gases. They are described by the limiting case of the equation of motion ([5]) in which 
one should replace the lengths A l and As by the single parameter R whose magnitude is of the order of the size of the 
vessel or the size of the cloud of the condensed atoms [1^. Of course, the term oc /i/ is excluded in this case because 
gauge field is absent, hence j> rda = 0. In this case, the torsion and writhe are conserved separately. 
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